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Abstract 

We construct a coarse-grained model of parallel actin bundles crosslinked by compact, globular 
bundling proteins, such as fascin and espin, necessary components of filapodial and mechanosensory 
bundles. Consistent with structural observations of bundles, we find that the optimal geometry 
for crosslinking is overtwisted, requiring a coherent structural change of the helical geometry of 
the filaments. We study the linker-dependent thermodynamic transition of bundled actin filaments 
from their native state to the overtwisted state and map out the "twist-state" phase diagram 
in terms of the availability as well as the flexibility of crosslinker proteins. We predict that the 
transition from the uncrosslinked to fully-crosslinked state is highly sensitive to linker flexibility: 
flexible crosslinking smoothly distorts the twist-state of bundled filaments, while rigidly crosslinked 
bundles undergo a phase transition, rapidly overtwisting filaments over a narrow range of free 
crosslinker concentrations. Additionally, we predict a rich spectrum of intermediate structures, 
composed of alternating domains of sparsely-bound (untwisted) and strongly-bound (overtwisted) 
filaments. This model reveals that subtle differences in crosslinking agents themselves modify not 
only the detailed structure of parallel actin bundles, but also the thermodynamic pathway by which 
they form. 
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I. INTRODUCTION 
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Parallel actin bundles are highly organized structures crucial to diverse range of cellular 
function, from mechanosensory specializations such as microvilli, stereocilia and neurosen- 
sory bristles to the highly dynamic filapodial protrusions of cell cytoskeletons 
assemblies share a common structural organization: axially-aligned actin filaments of uni- 
form polarity, densely-assembled into an ordered hexagonal array and interspersed with a 
crosslinking array of actin bundling proteins. Multiple bundling proteins have been identi- 
fied from parallel actin bundles in vivo, though the type and composition bundling proteins 
varies significantly between different cell types 4|]. Primary examples of actin bundling 
proteins, fascin and espin, are known to be integral components of filapodia and stere- 
ocilia bundles Qj, respectively. It is believed that the array of multiple bundlingproteins 
affords cells the ability to form actin bundles with variable properties, such as size |8L M and 
rigidity [1CH121]. though little is understood about how distinct features of bundling proteins 
specifically modify the assembly of actin filaments into bundles. 

Structural studies of in vitro bundles 



18| suggest that a key aspect of the forma- 



tion of parallel bundles is the ability of crosslinking proteins to modify the twist of actin 
filaments. In bundles, the helical symmetry of constituent filaments is modified from its 
native -13/6 geometry: a left-handed helix rotating through 6 turns per 13 monomer repeat. 
This native geometry is poorly suited for the hexagonal symmetry of the array, which favors 



co-registry of crosslinked monomers on neighbor filaments 15|. Electron diffraction studies 



of fascin-crosslinked bundles reveal that filaments are overtwisted to a —28/13 symmetry, 
corresponding to a change of roughly - 0.01 monomers /turn, a distortion which is consistent 

tin 

with more recent observations of espin- mediated bundles [13, H8|. Despite the apparently 
similar structural change induced by crosslinking, fascin- and espin-mediated bundles ex- 
hibit a dramatically different sensitivity to concentration of available crosslinkers in in vitro 
systems. Based on small-angle x-ray studies, Claessens et al. found that overtwist of fila- 
ments in fascin-mediated bundles is sensitive to the concentration of available crosslinker, 
with helical filament symmetry varying continuously from native to fully overtwisted sym- 
metry In comparison, a recent study of espin cross-linked bundles found that above a 
critical concentration of crosslinker, bundles lock into the fully overtwisted state, with little 
or no further sensitivity to espin concentration [jij. Evidently these two compact, globu- 



lar bundling proteins, fascin and espin, are capable of forming bundles of the apparently 
same overtwisted structure, though in each case the fully-bundled state is approached via a 
different pathway of states at intermediate crosslinker concentrations. 

This rich phenomenology raises a number of questions about the role of filament twist in 
the assembly mechanism of protein-mediated bundles. What is the mechanical cost associ- 
ated with distorting bundled filaments from their native geometry, and how does this cost 
effect the bundling transition? What is the role of the torsional rigidity of the crosslinking 
bonds themselves? Most puzzling, what is the nature of the states of intermediate twist 
observed for fascin-mediated bundles? 

In this article, we analyze a quantitative model that describes the complex interplay be- 
tween the optimal geometry required by crosslinking in actin in hexagonal bundles and the 



cost of distortin, 
posed in ref 



g filaments from their ideal helical symmetry. Based on a lattice model pro- 
we identify a unique crosslinked-bundle groundstate with -28/13 symmetry 
that allows for an optimal number of "ideally" oriented crosslinking bonds. We study the 
thermodynamic transition from untwisted, unbound filaments to fully-bound, overtwisted 
bundles driven by increasing the concentration, or chemical potential, of available crosslink- 
ers. A coarse-grained model of parallel bundles, maps the linker-induced overtwist of actin 
filaments onto a commensurate-incommensurate phase transition. We find that this transi- 
tion takes place by a surprisingly complex coherent restructuring of filaments in the bundle, 
in which bundles possess localized bands, or domains, of native (-13/6) and overtwisted (- 
28/13) filaments. The overtwisting of bundled filaments then proceeds as the fraction of 
overtwisted bundles increases continuously from 0, in the absence of crosslinkers, to 1 in 
excess of available crosslinkers. Owing to the fundamental role of elastic distortions in this 
model, the bundling transition is found to be extremely sensitive to the flexibility of the 
crosslinking bounds. For sufficiently rigid crosslinkers, bundles pass to the overtwisted state 
via a sharp, thermodynamic transition; while for bundles held together by relatively flexible 
linkers, a smooth transition to a maximum state of twist is predicted. The primary conclu- 
sion of this study is that differences in the sensitivity of the bundling transition to crosslinker 
concentration observed from fascin- and espin-mediated bundles derive from distinctions of 
the flexibility of crosslinker to actin bonds. Hence, we have identified the flexibility of the 
crosslinking bonds provided by bundling proteins as a key parameter controlling not only 
the structure, but the process, by which parallel actin bundles are formed in different cell 




FIG. 1. A schematic picture of two actin filaments crosslinked by a bundling proteins (purple) 
and its lattice model representation in a box, where G-actin monomers are described by a set of 
vectors, S^. The top view of layer i is also shown with the angular deviation 5(j)n of the monomer 
from the lattice direction Dy. 

types. 



II. LATTICE MODEL OF CROSSLINKING IN PARALLEL ACTIN BUNDLES 



We model a bundle as a parallel, hexagonally-ordered array of actin filaments with fixed 
nearest neighbor spacing D f=s 17 nm, consistent with structural observations 

The 

positions of the monomeric, G-actin are described by a set of vectors, aSi/, that point from 
the center line of the ith actin filament of the lattice to the center of the £th monomer 
along the filament, where a ~ 3.75 nm roughly the diameter of G-actin (see Fig. [TJ). In 
the native twist state, these vectors precess around the centerline of filaments at a constant 



angular rate of Uq = 127r/13 per monomer (i.e. 6 rotations per 13 monomer repeat) 19]. 



We describe torsional distortions with the following elastic energy 
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where C is the torsional elastic modulus of actin filaments [20[, and A<^ = <fti/+i — <ftij is 
the rotation angle between successive monomers along the ith filament, with the angle 
of the £ih monomer direction in the plane of lattice order. 

In our model, crosslinking between neighboring filaments in the bundle occurs between 
pairs of monomers at the same vertical layer, £, shown schematically in Fig. Q] . Because 



globular bundling proteins like epsin and fascin have compact size, ~ 5 — 7 nm in diameter 



21], in comparison the lattice between filament in the bundle ~ 15 — 20 nm, crosslinking 



occurs preferentially when monomers on adjacent filaments are minimally separated from 
one another. Yet, due to incommensurate helical symmetry of actin filaments, bundles are 
forced to accomodate crosslinks with some degree of "misfit" between crosslinked monomers. 
We describe this effect with the following simple elastic model for enthalpy of crosslinking 
between monomers, 

Eunding = ^ n ijti - € b + ^(5(j)\ t + 5^) , (2) 

where the sum is over neighbor filaments, (ij). In this formula, is 1 if the bond between 
i and j at £ is occupied, and if empty, e b describes the energy gain of a perfectly-aligned, 
crosslinking bond between 2 monomers, and the final term in the parentheses describes 
the energetic cost of distorting the bond from its ideal geometry. Here, 54>i t e is the angular 
deviation between Si/ and D^, the lattice vector separating % and j. Finally, U is a parameter 
describing the "elastic cost" of distorting the aligned bond between monomers. For example, 
if this cost could be described purely in terms of a simple linear spring energy, fc^A^ — 
Ao) 2 /2, which penalizes changes in length, A^, of the monomer-monomer separation from 
a zero stretch length, Ao = D — 2a the equilibrium size of crosslinks, the effective elastic 
parameter in ([2]) becomes U ~ 2/^a 2 : . 

The bundling of actin is sensitive to the concentration of free crosslinkers in solution. 
We therefore study the thermodynamics of crosslinking at a fixed chemical potential, /i. 
This accounts for the equilibrium free energy cost of removing a free crosslinking protein 
from solution and adding it to a bundle, and therefore, \i is related to the free crosslinker 
concentration by c/ ree oc e M//fesT . 



III. IDEAL CROSSLINKING GEOMETRY IN BUNDLES 

To describe the overtwist transition of parallel bundles, it is necessary to understand 
the optimal geometry of highly crosslinked bundles, as well as the low energy pathways to 
this state from the native actin geometry. First, we briefly demonstrate the structure of 
optimally-packed actin bundles in terms of geometric considerations (see supporting ma- 
terials for full details). The model described above is highly frustrated, a generic feature 
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of hexagonally-organized filament assemblies, well-studied in the context of counterion me 



diated biopolymer bundles 
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a-|24j as well as helically-ordered phases of DNA 25N27|. 
Here, we consider the configurations for which crosslinking bonds are perfectly aligned to 
the lattice directions (that is, = for all = 1 ) and for which configurations the 
number of perfectly oriented bonds is maximal. 

Perfectly aligned configurations require a subset of the actin monomers to align with 
a sixfold lattice direction. We construct structures which have alternating sequences of 
sections of —6m/n symmetry - each successive monomer is rotated by —27m/ (6m). Here, n 
and m are integers so that mth monomer lines up with lattice direction of 27m/6, allowing 
for a perfectly aligned crosslinking bond to form. It is not difficult to show that among 
these commensurate helical geometries the —24/11 and —30/14 structures are particularly 
close the native geometry of actin, differing only by 0.69% and 1.1%, respectively, in terms 
rotation angle per monomer. This proximity to the native geometries confers upon them an 
especially low twist cost among all possible ideal crosslinking states. To determine the bundle 
structure with the maximum number of bonds, we therefore considered periodic states with 
a composite symmetry, possessing N4 numbers of 4- monomer sections with -24/11 symmetry 
and N 5 numbers of 5-monomer sections with -30/14 symmetry. To construct crosslinks, it 
is not sufficient to consider aligned monomers to lattice directions from a single filament, 
as crosslinking requires the co-orientation of monomers on neighbor filaments at the same 
vertical layer along the filaments. Hence, it is necessary to consider the three-dimensional 
geometry of possible multi-filament structures arrayed on the hexagonal lattice. 

Based on an extensive numerical search of composite —24/11 and —30/14 structures up 
to 102 monomers per repeat length, we find a maximum crosslinking density for N4 = 2 
and N$ = 4, which has 6 crosslinks along every 28 monomer length of actin filament (see 
Fig. E]). Notice that this composite structure has -4iV 4 (11/24) - 5iV 5 (14/30) = -13 net 
turns per 28-monomer repeat. The overall symmetry and bond/monomer stoichometry of 
this ideal geometry are in perfect agreement with careful structural studies of overtwist actin 
bundles, formed by both espin and fascin crosslinkers, which also have an overtwisted —28/13 



structure and crosslinks spaced at 4- and 5-monomer intervals along filaments [151 ] . Though it 
has non-hexagonal symmetry, the composite filament-bond structure of this perfect packing 
geometry can be repeated to construct a parallel bundle of arbitrary size. This unique 
crosslinker geometry serves as the overtwist groundstate of our model. 
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FIG. 2. The unit cell of maximally crosslinked-bundle groundstate in the plane of hexagonal order. 
The repeat unit contains a single filament geometry translated and rotated vertically to obtain the 
geometry of the 4 filaments. Numbers represent the vertical layers of co-oriented monomers with 
the lattice direction, and gray bars, crosslinkers. On the right, the 28-monomer vertical repeat 
geometry of 2 neighbor filaments, highlighting definition of the coherent rotation of filaments 
from the ideal crosslinking geometry. 



IV. UNTWISTING OVERTWISTED BUNDLES: COARSE-GRAINED THEORY 



Having identified the limiting geometry of unbound filaments (13-fold helical symmetry) 
and fully bound filaments (28-fold helical symmetry), we consider the thermodynamic pro- 
gression of filament twist as crosslinker density in bundles increases. Although the detailed 
structure of the overtwisted -28/13 structure is somewhat complex, the underlying screw- 
symmetry of actin filaments and the in-plane periodic order imbue this state with a rather 
simple symmetry under coherent rotations of each filament by 2n/28 around its axis. A 
27r/28 rotation of each filament in the ideal crosslinking configuration in Fig. [2] followed by 
a rearrangement of monomers and bonds within the unit cell recovers an equivalently ideal 
geometry, with 6 perfectly aligned crosslinks per 28 monomers (see supporting information). 
Hence, the ultimate function of the complex pattern of crosslinking bonds is to lock the bun- 
dle into a -28/13 twist-symmetry and constrain the azimuthal orientation of this structure 
to within one of 28-fold bonding free energy minima. Competing with this tendency is the 
intrinsic torsional elastic energy of filaments which favors unwinding the overtwisted state 
to the -13/6 symmetry, making it costly for the bundle to maintain a 28-fold commensurate 
bond geometry along its length. 

7 



To model the free energy gain of crosslinking, we construct actin bundles in the fully- 
overtwisted state with a net degree of twist, Qq = 27r(13/28), and consider the low-energy 
distortions that untwist the 28-fold commensurate geometry of bundles as crosslinkers un- 
bind from bundles. The analysis is based on a coarse-graining of the model described by 
eqs. ([T]) and ($Z§. In particular, we decompose the monomer orientations in terms of two 
angular deviations from homogenuously overtwisted filaments: (j)^, which describes short- 
lengthscale monomer relaxations within a 28- monomer repeat length, and describing the 
slow, coherent rotations of filaments on much longer length scales. In terms of the angle a 
monomer direction makes in the plane of hexagonal order, we define, 



where we restrict the lengthscale reorganization to sum to zero rotation within a 28-monomer 
repeat length, Yfe^e^ = so that the net rotation of filaments away from the overtwisted 
state is $£. As a description of the long- lengthscale structure of bundles, <3>£ serves as the 
order parameter, fully describing the underlying state of our model: overtwisted states com- 
mensurate with the ideal crosslinking geometry of the bundle correspond to = 27rm/28 
for any integer, m. 

Due to the separation of length scales between the deformations described by and <pi,i 
the elastic twist energy approximately decouples these degrees of freedom, 



where 5uq = u>o — Qq = 27r/364 is the overtwist distortion per monomer from the native 
to the homogeneously-twisted -28/13 state and we have implicitly assumed that is 
approximately constant over the length of a 28-monomer repeat. 

The final step of our coarse-graining minimizes the twist and binding free energy, eqs. 
(jH) and (j2|), over the distributions of bonds, n^, and 0^, the short length-scale angular 
reorganizations for a given value of $. We perform this minimization by analyzing a 28- 
monomer repeat of the ideal bonding configuration Fig. [5J requiring a net rotation of this 
structure by $, and finding the subset of the 6 perfectly aligned crosslinks for which net 
free energy of binding, including the cost of elastic distortion, is minimal (see Appendix IA1) . 
The result is the binding free energy per monomer, V($), that depends only on mean value 



>i,e = 



n £ + $ t + <j> ift , 



(3) 




(4) 
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of $ within a 28-monomer repeat (the coarse-grained unit of our model), 

V@) = min nb , m [ - y (/i - + ^ ($ + ^) *] /28, (5) 

where n& is the number of crosslinked monomers per 28-monomer repeat, which varies from 
n& — (unbound) to = 6 (fully bound) and the minimization of m reflects rotational 
symmetry of the binding free energy, V($ + 27r/28) = V(&). Here, /i c = — e& + S/i nb , 
where 5fi nb represents the excess torsional elastic energy per bond needed to distort the 
homogeneously overtwisted filament into a state where the n& bonds are perfectly aligned 
with bond directions. This torsional cost represents an offset to the binding free energy of 
the state with n& monomers proportional to C that increases with number of bonds: 5^o = 0; 
Sfji! = 0; 6 & = 0; Sfi 3 = 0.00058C; fyx 4 = 0.00087C; 5/x 5 = 0.00191C; and 5fi 6 = 0.00261C. 
The term proportional to $ 2 represents the resistance of the structure to rotations from 
ideal crosslinking geometry, 

_ n b U 

This elastic response of the crosslinking array to coherent rotation is straightforward to 
understand in the small and large U limits. The rotation of filaments from the ideal bond- 
ing state requires either the bond orientation - as parameterized by S<fri/ for crosslinked 
monomers - or the torsional state of the filament within the 28-monomer repeat to adjust. 
When (7/C < 1 and linkers are more flexible than the filaments, this torsional load is carried 
by the flexibility of the crosslinks themselves, hence, K nb oc n^U. For very rigid linkers, 
!7/C > 1, bound monomers are pinned to the lattice directions the rotation of the filament 
section is accomplished instead by a twist distortion of the monomer segments neighboring 
bonds, so that K rib oc n^C. 

Shown for flexible and stiff crosslinkers in Fig. El V($) functions as the "rotational 
potential" describing the free energy preference for the filament structure to lock-in to a set 
of preferred torsional configurations due the favorable number and arrangement of bonds 
in overtwisted bundles. By definition the coherent rotation varies slowly on the scale of 
monomers. We take the continuum limit of our model, $(^) = and write the final form 
of the free energy of a bundle with rif filaments as, 

-L 



bundle 




FIG. 3. The free energy gain per monomer of crosslinking as a function of <£, the angle of coherent 
rotation of filaments in the bundle for different values of the linker chemical potential, eq. ([5]). 
Highlighted in the for the rigid linker case, U = 3C, are regions of the rotational potential where 
a 28-monomer repeat maintains 0, 2, 4 and 6 bonds. The maximal number of bonds occurs for 
commensurate geometries with <I> = 27rm/28. 



Written as such, the coarse-grained free energy highlights the essential frustration of parallel 
actin bundles. The first term in the integrand is minimized when d^/dz = 5u>o and the 
filaments in the bundle revert to the native, 13-fold helical geometry. Competing with this 
is the rotational potential, which is minimized by a constant value $ = 27rm/28. The relative 
importance of these competing effects is sensitive to fi, which largely dictates the depth of 
V($), but also U and C which together determine the relative stiffness of filaments and the 
pinning of rotational potential. 

The effective model for parallel actin bundles, eq. (J7D, is known in condensed matter 
contexts as the Frenkel-Kontorowa model, employed in the study of incommensurate, one- 
dimensional solids 30|. The structure and thermodynamics of the minimal energy ground 
states show a complex dependence on the degree of incommensurability, 5uq, and relative 
strength of the potential pinning the solid in the commensurate state (here the overtwisted - 

10 



28/13 bundle) to the elastic energy of the incommensurate state (the native -13/6 symmetry). 



The free-energy minimizing solutions are described the following differential equation 

C fd§\ 2 



29|, 



2 UJ =n*)-no) + e, (8) 

where e is a non-negative parameter that specifies the entire rotational structure, along 
the bundle. Minimizing the bundle free energy (j7|) of this class of solutions yields an equation 
for e which corresponds to the mean structure of the bundle, 

5u = / d$y/V($) - V(0) + e . (9) 



'2Ctx Jo 

Even in the absence of thermal fluctuations, this model has a complex dependence on the 
binding free energy. For sufficiently strong pinning potentials, the lowest energy solution 
becomes $ = and e = 0, indicating that the bundle has locked- into the commensurate 
phase, here -28/13 overtwisted structure. Below a critical depth of the pinning potential, 
solutions with e > exist indicating that $(£) has an inhomogeneous solution, which grad- 
ually unwinds to the native state. This corresponds to the incommensurate phase of the 
Frenkel-Kontorowa model. As a measure of the average rate of rotation of filament struc- 
ture, we define the length £ as the length along which the filament geometry unwinds by 
27r/28, from one minimum of V(<&) to the next. From eq. (jHJ) this length, measured in 
monomer number, is computed from the integral, 

fc r 2w/28 d$ 

C=\- . (10) 

V 2 J ^($)-y(0) + e 

This length is related to the mean rate of filament twist by, 

(A0> = Q - —C-\ (11) 

In the following section, we analyze the behavior of this order parameter, as well as the 
detailed structure of parallel bundles in terms of the inhomogeneous solutions for filament 
rotation, $(f). 

V. OVERTWIST TRANSITION 

The overtwist thermodynamics of crosslinked parallel-actin bundles predicted from the 
coarse-grained model is shown in Fig. HI which shows the mean filament twist in terms of 
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FIG. 4. The twist-state phase diagram - contours of ( A0) - in the parameter space of (fi + £fe)/C 
and J7/C. The evolution of the stable twist state from -28/13 (overtwist) to -13/6 (native) state is 
observed. 

crosslinker chemical potential, /x, and the stiffness of crosslinking bonds, U. For \x < — e&, 
crosslinks are not favored in the bundle, and hence V($) = 0, indicating no thermodynamic 
preference for twist, hence, in this region (A0) = oj , the native state of actin geometry. 
For all values of fi > — e&, crosslinking is favorable and the incommensurate geometry of the 
-28/13 ideal crosslinker favors an overtwisted actin geometry with (A0) > uq. 




FIG. 5. The mean rate of twist (A0) per monomer, from ujq to fio , is plotted as a function of the 
chemical potential (// + £fe)/C for the various values of J7/C. The dotted line is for critical value 
of (fi + e&) = 0, where the minimum in the binding free energy potential abruptly changes from 
V(0) = to V(0) < 0. 
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The sensitivity of filament twist to fi + e& is dramatically different for stiff and flexible 
crosslinking bonds. Shown in Fig. [5] is the //-dependence of filament twist plotted for linker 
stiffnesses ranging from more rigid (U > C) to more flexible (U < C) than effective torsional 
resistance of 2 successive monomers. For rigid linkers, we predict that filament twist initially 
increases continuously with increased /i near the onset of crosslinker binding /i > — e&, but 
upon approaching a critical value of the chemical potential, nch very rapidly overtwists to 
the -28/13 structure. For \i > /ici, the filament geometry locks into (A</>) = and exhibits 
no further sensitivity to the availability of crosslinking proteins. 

This singular dependence of (A0) on /i+ej, is the signature of commensurate-incommensurate 
(CI) phase transition of eq. (J7|), marked by a divergence of the distance over which $ rotates 
by 27r/28, £ ~ — ln(//ci — /f)- The extreme sensitivity of bundle structure to crosslinker 
chemical potential near the CI phase transition derives from the highly cooperative change 
of symmetry of low energy state of rigidly crosslinked bundles mediated by the elasticity 
of the filaments and the array of crosslinking bonds in the bundles. The value of /ici is 
determined by the /x at which e — > from eq. ([9]), is shown in Fig. HI In the limit of infinite 
linker rigidity, nci(U 3> C) — > — e& + 0.0038C, and this critical value shifts to larger /i as 
U is decreased. In contrast to the stiff linker limit, below a critical value of linker stiffness, 
[/* = 0.512C, the overtwist shows no CI transition, and (A</>) shows a considerably reduced 
sensitivity to /i. 

What accounts for the distinction between bundling by rigid and flexible crosslinkers? 
In our coarse-grained model, these differences in thermodynamic behavior ultimately derive 
from distinct features rotational free energy potential, V (<§>), describing the sensitivity to 
the free energy gained by crosslinking to rotational state of actin filaments in the bundle. 
The thermodynamic preference to lock into the -28/23 structure can be crudely under- 
stood in terms of relative cost overtwisting of the filament to a constant $ = state, 
C(Sujq) 2 /2 per monomer, and free energy cost of rotating the ideal filament structure from 
the commensurate state, roughly corresponding to the depth of the rotational potential, 
AV = V(ir/28) — V(0), describing the thermodynamic cost of the non-ideal crosslinking 
geometry. From eqs. (EJ) and (EJ), the narrowness of the minima of V(<&) is determined by 
K nb , which increases monotonically with increased linker stiffness. As shown in Fig. [31 rigid 
linkers only allow favorable crosslinking for a narrow range of rotations from the commensu- 
rate geometry, and hence, rotating away from $ between minima necessarily forces filament 
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FIG. 6. The mean number of crosslinkers, {rib), P er 28-monomer repeat unit. The solid curve indi- 
cates the commensurate-incommensurate transition, while the vertical lines at fj, + €b = 0, 0.0017C 
and 0.0061C, denote abrupt transitions in minimal-energy number of crosslinkers as described by 
eq. ©. 

to release its crosslinks, AV ~ (3/28) (/i + e 6 — <5/i 6 ). Thus, for sufficiently large \i this penalty 
outweighs the cost of overtwist, and rigid crosslinkers always lead to bundles locked-into the 
commensurate state. For flexible linkers, V($) is a shallow function of rotation angle, and 
all rotation angles up to $ = n/28 allow filament segments to maintain 6 crosslinks, albeit 
somewhat stretched from the ideal geometry. In this case, AV m (3/28)£/(7r/28) 2 , which 
is determined entirely by the elastic cost of stretching linkers as no crosslinks need unbind 
to overcome the free energy barrier. Thus, sufficiently flexible linkers allow the bundles 
to accommodate a maximal number of favorable crosslinking proteins even in the absence 
of a dramatically overtwisted structure. Therefore, we see that {rib) the mean number of 
crosslinkers per 28 monomer repeat (see Appendix [B]), shown in Fig. is predicted to 
increase to 6 crosslinks for large e&, independent of crosslink flexibility. 

VI. STATES OF INTERMEDIATE TWIST 



Upon crosslinking, a parallel actin bundle undergoes a complex structural change, from 
a state of filaments possessing a native helical geometry, to a fully-bundled state with an 
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ultimate geometry that is sensitive to the flexibility of the crosslinking bonds themselves. 
The one-dimensional, coarse-grained model of eq. (jTJ) predicts that a bundle progresses 
between these limiting structures via a rich pathway of intermediate states of inhomogeneous 
filament twist. Shown in Fig. [7J are minimal energy solutions of eqs. OH]) and (JH]) for the 
lengthscale, coherent rotations of bundles filaments, <&(£), for U = 2C. 

Beginning in the commensurate state in excess of available crosslinkers (/i > /xci)> the 
bundle filaments lock into the ground state of binding free energy density, $(^) = 0. For fi 
just below this critical value, a weaker rotational potential does not hold the filament in the 
commensurate overtwist state along the entire filament length. Instead, filaments untwist by 
rotating from one commensurate orientation to the next, say from 27rm/28 to 2ir(m + l)/28, 
by way of rapid "jumps" in In the language of incommensurate solids, these anglar 

jumps are known as solitons or discommensurations. In this model a discommensuration 
spans the cross section of the bundle, representing a domain of nearly native filament geom- 
etry, d(p/d£ ~ 5uo- Over a relatively short span of roughly ~ 15 — 20 monomers $ rotates 
between two nearby minima in the binding free energy. For rigidly crosslinked bundles, the 
mean number of crosslinks in these domains is significantly reduced from n& = 6 per 28 
monomer sections, as crosslinks are predicted to unbind for at the free energy maximum of 
V($) (see Fig. Ej). 

The inhomogeneous twist and crosslinking structure of parallel bundles is shown skemat- 
ically in Fig. The minimal free energy bundle configurations are described by periodic 
solutions for in which discommensurations have an equilibrium spacing, C, along the 

bundle. As this length becomes very large close to the /jLqi, the regions between discommensu- 
rations represent sections of overtwisted, maximally crosslinked bundles. Hence, a surprising 
prediction of this coarse-grained model is that states of intermediate twist are constructed of 
alternating domains of sparsely-bound, native filament geometry and strongly-bound, over- 
twisted filament geometry. The net rate of filament twist, (A0), is predicted to increase as 
the relative proportion of overtwist to native twist filament geometry increases, as further 
portions of bundle are converted to highly-crosslinked, -28/13 geometry. As shown in Fig. [JJ 
when fi is reduced well below /ici, C decreases until successive discommensurations merge, 
relaxing the filament to its homogeneous state of native geometry. 
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FIG. 7. The profile of coherent rotation $ with respect to the I monomer is plotted at the various 
values of fx + e^. On the right, the inhomogeneous structure of parallel bundles, composed by 
alternative domains of sparsely-bound native (blue) and densely-bound, overtwist (red) symmetry, 
is shown schematically. 

VII. DISCUSSION 

When actin filaments are bundled in parallel arrays by action of compact, globular 
bundling proteins, the process of crosslinking affects complex structural change of the fil- 
aments themselves, an overtwisting from a native -13/6 helix to -28/13 symmetry. This 
structural change is required to maximize the number of co-oriented monomers on neighbor- 
ing filaments, themselves arrayed on a hexagonal lattice. The necessity to overtwist actin 
filaments upon binding of crosslinking proteins leads to a complex thermodynamic depen- 
dence of bundle properties - mean twist and bound crosslinker density - on the availability 
of crosslinkers. As it is the bonds themselves that mechanically distort the helical structure 
of bundled filaments, we find that the thermodynamics of the bundling process is extremely 
sensitive to the rigidity of the crosslinking bonds. A mapping of a coarse-grained free energy 
of bundles onto the Frenkel-Kontorowa model of incommensurate solids predicts that rigid 
crosslinkers (U > £/*) rapidly overtwist filaments via a second-order phase transition, that 
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takes place over a very narrow range of crosslinker chemical potential. In contrast, flexible 
crosslinks are predicted to be much less efficient in providing the necessary overtwist of 
filaments, leading to a bundling transition that takes place over a much broader range of 
chemical potential, and consequently, over broad range of free crosslinker concentrations. 

This distinction in crosslinker behavior between rigid and flexible crosslinkers predicted 
here is entirely consistent with the observed differences between crosslinking behavior be- 
tween the globular crosslinking proteins, espin 17, 18] and fascin jsl, [isj]. Based on our 



present study, we attribute the respectively abrupt and continuous dependence of filament 
overtwist observed for espin- and fascin-mediated bundles, not to differences in linker affin- 
ity, but to differences in flexibility of the crosslinking bond. This suggests that, in part, 
different cells form bundles by way of different bundling proteins, in order to regulate not 
just the structure of bundles, but also self-assembly pathway of bundles and consequent 
sensitivity of these bundles to modulations in the availability of bundling proteins. By al- 
tering the flexibility of the crosslinking agents alone, the formation of bundles can be tuned 
from a highly-cooperative, switch-like dependence to non-cooperative, continually-varying 
dependence on crosslinker concentration. 

A critical value of crosslink stiffness, C7* = 0.512C, separates the stiff linker behavior, 
which exhibits a CI transition, from the flexible linker behavior, exhibiting no overtwist 
phase transition. The twist modulus of actin is measured to be C ~ 2.8 x 10~ 17 J ~ 6900 k^T 



(units of angular distortion per successive monomer pair) 20] . If we attribute the elastic 
cost of distorting binding directions to a change of crosslinker length, we may estimate 
an effective spring constant for linkers with this critical stiffness k* ~ 500 pN/nm, based 



a ~ 3.75 nm, half the diameter of actin 19]. This is considerably more rigid than the bonds 



provided by the relatively extended Crosslin 
larger bundling proteins, filamin and a-actin 



ters in acto- myosin bundles [31], or the much 



32] . all of order 0.1 — 10 pN/nm. The stiffness of 



crosslinking bonds provided by small, globular bundling proteins is not a well characterized 
quantity, though it seems reasonable that U ~ C, since the torsional flexibilty of the actin 
filaments requires distortion of specific protein-protein bonds between actin monomers that 
may not be wholly unlike the distortion of the bonds between actin monomers and bundling 



proteins of roughly the same size. Based on a mechanical model of crosslinker shear 33] . 
measurement of the effective ben ding stiffness of parallel bundles have been used to infer a 
shear stiffness of fascin crosslinks 10j that is 3-4 orders of magnitude below this estimate for 
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fc*, suggesting that fascin should fall well within the flexible linker regime. Notwithstanding 
the poor understanding of the mechanics of crosslinking at the protein scale, the classification 
of fascin as a flexible bundling protein is consistent with its role as a primary crosslinker 
in filapodial bundles Q. Fascin has been observed to transiently bind and unbind along 
filapodia in vivo, diffusing at a remarkable rate within bundles |34j. In the context of our 
present study, we note that flexible linkers bind and unbind in an essentially independent, 
non-cooperative manner. In contrast, we expect the kinetics of rigid crosslinkers to be 
dramatically slower than flexible linkers due to the cooperative organization of many linkers 
and coherent restructuring of filaments required for each additional rigid crosslinking bond. 

We conclude by discussing the role of two effects not included in our model crosslinked 
parallel actin bundles: thermal fluctuations of monomers and crosslinks and the global twist- 
ing of finite bundles. In the present model, the complex thermodynamic binding properties 
of crosslinking proteins derives purely from an elastic frustration between optimal linker 
and filament geometries, neglecting thermal fluctuations of monomers and bonds. In ref. 



181 ] we studied this lattice model of parallel bundles in the presence of strong thermal 
fluctuations of filament orientation and crosslinker position. In this regime, the thermo- 
dynamic distinction between by rigid and flexible bundling transition is maintained: rigid 
crosslinkers affect a phase transition between native and overtwist, while the flexible linkers 
continuously overtwist filaments upon increasing crosslinker binding. Given the high tor- 



sional modulus |20| , one expects a torsional persistence length of free, unbundled filaments, 
ft = 2C/k B T ~ 13,000 monomers defined by (cos(<^ - <pse+e )) = cos(w ^)e- |w|/5 *. On 
the scale of the coarse-grained repeat lengths of our model, 28 monomers, thermal fluctu- 
ations of filament twist are therefore extremely modest. However, we do expect that the 
entropy associated with distributing crosslinkers at different positions in bundles to quan- 
titatively modify the predictions of our present study, although modestly. In particular, 
the non-analytic dependence of the binding free energy, V(<&), on chemical potential and 
filament rotation is smoothed out by the equilibrium distribution of crosslinks among the 
many competing crosslinking geometries (see supporting information). Therefore, certain 
sharp features of binding thermodynamics predicted by this "zero temperature" theory are 
somewhat smoothed out when considering these fluctuations. In particular, the number of 
bonds in the minimal free energy states changes abruptly from n& = to 2 at \i + €& = 0, 2 
to 4 at fi + e b = 0.0017C and 4 to 6 at /i + e& = 0.0061C, leading to kinks in the predicted 
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equations of state for (rib) and (A0) vs. \i (Figs. [5] and [6]), which only appear in this limit. 
Nonetheless, we find that presence of these fluctuations does not eliminate the sharp CI na- 
tive/overtwist transition that occurs for linkers above a critical stiffness, £7* = 0.512C, where 
d(<p)/dfi diverges. That is, the sharp/smooth overtwist transition for rigid/flexible crosslink- 
ers is robust feature of the geometrical frustration in parallel actin bundles, insensitive to 
the presence of thermal fluctuations of bundles. 



In our model, we considered parallel actin bundles of essentially unlimited length and 
width, focussing on the internal reorganization of filaments within the bundle. Cells maintain 
a careful control over the size of bundles, most notably in the mechanosensory hair bundles 
of the cochlea 35|. Understanding the physical mechanisms that underlie the in vivo control 



of bundles on multi-filament lengthscales (~ 0.1 — 10 /im) remains an outstanding puzzle. In 
vitro studies of bundle formation observe that the lateral diameter is indeed sensitive to the 
concentration of crosslinkers in fascin/actin solutions 



Claessens et al. suggest that the 



observed coincident overtwisting of filaments and growth of bundle diameter upon increased 



fascin binding implies that overtwist plays a role in 



actin bundles 
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imiting the the lateral assembly of 



40] demonstrate that a tendency of 



9j. Indeed several theoretical studies 
filaments to globally twist around the central axis of the bundle leads to thermodynamic 
frustration that ultimately limits the equilibrium diameter of bundles. This mechanism is 



41} , which are clearly observed 



believed to play a role in the self-limited size of fibrin bundles 
to globally twist in electron microscopy studies. A similar mechanism may be at work in 
fascin-mediated bundles, as the intra-filament and inter-filament twist of helical filaments 
are geometrically coupled j^]. To understand the relationship between crosslinker binding 
and bundle size, it is therefore necessary to consider a more complex model of protein- 
mediate parallel bundle formation, in which filaments trade the elastic cost of filament 
overtwist for the elastic cost of globally twisting the entire bundle. Though it is plausible 
that the tendency to overtwist individual filaments also implies a tendancy to twist the 
entire filament lattice in a bundle, it remains to open to question whether the elasticity of 
crosslinking bonds is of sufficient rigidity to mechanically bend filaments into the complex 
superhelical structures that frustrate and limit self-assembly. 
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Appendix A: Rotational potential 

We drive the explicit form of the potential V($) in eq (JSJ) by considering a 28-monomer 
repeat length of filament in which rib of the ideal bond structure form with neighboring 
filaments. We describe the angles of this initial state with n& monomers aligned to six- 
fold lattice directions by the angles <p™ b , and each of these states has a composite -28/13 
symmetry. As bonds are elastic and commensurate orientation varies along the bundle 
length, the perfect orientation of these bonds will relax to a low energy configuration. 

To determine the relaxation of orientation within a 28-mononer unit cell, we introduce 
two angular deviations: 5<frb, the crosslinker deformation from the perfectly aligned state at 
the crosslinked layers and S(pf the adjustment of unbound, or free, monomers. At layers 
bound by crosslinks, say £q, the free energy of bonding per monomer from eq. (J2J) is simply 
— (/i + e&)/2 + [7(50b) 2 /2, since each crosslink is shared by 2 filaments. Additionally, the 
twist elastic distortion between the £q and £q ± 1 monomers from eq. ([T]) contributes, 
C(A$£ + 5(f) b - 6(f) f ) 2 /2 + C(A0" o 6 +1 - 5(j) b + 5(f) f f/2 to the free energy of binding, where 
0" b = <f)™ b — Qq£. The twist elastic energy of the remaining unbound monomer sections 
is simply, C(A<j)™ b ) 2 /2. Summing these contributions, we have the free energy of binding 
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within a 28-monomer section per filament in the n b bond state, 

F%(6<j», 50/) = y (t/(50 fc ) 2 + 2C(50/ - 50 b ) 2 - (,x + e 6 )) + ^ f>0^) 2 . (Al) 

Since $ is the coarse-grained coherent rotation of all these deformation in a 28-monomer 

unit, we have the following relation, 

= n b S(f) b + (28 - n b )5<f> f 

28 1 J 

Minimizing eq. (lAlj) with respect to the monomer adjustment angle, 50b, for fixed $ we find 

the minimal energy bond angle deviation is 

$ 

(A3) 



u.(l _ m>)2 — I — x 

2C \ 28^ ~ x 



Note that in the limit of rigid crosslinkers, when U ^> C, the bound monomers are pinned to 
the lattice directions, 50;, = 0. The resulting minimal free energy of the n^-monomer state 
rotated to an angle, $, is 

^(*) = ^ 2 -y(/^-^K)) (A4) 

where K nb is defined as in eq. (jSJ) and ii c {n b ) = — e^ + C ^ 2 f :O (A0" i, ) 2 /n b . The rotational po- 
tential is determined by minimizing F 2 "g ($) over all possible bond configurations within the 
28 monomer repeat for a given value of $, according to eq. (J5J). Though are a total of states 
considered within this class, for a given bond number a single bond configuration minimizes 
the elastic "offset" energy, C Y^oi^^T) 2 ( see table II in the supporting materials). 



Appendix B: Mean crosslink number 

At a given coherent rotation $, the optimal number state of crosslinkers, n b ($>), is de- 
fined as the number state which minimizes eq. ( 1A4|) over all the states of perfectly aligned 
crosslinking, from n b = to n b = 6. As shown in examples of stiff crosslinkers in Fig. [3j the 
minimized potential V(<&) possesses the subset of certain number states over the periodic 
range of $. 

In the minimal free energy configurations, $(£) rotates through multiple angles along the 
bundle corresponding to the solution of eq. (jHJ), 

i =L • (B1) 
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We calculate the mean number of crosslinkers in a 28 monomer repeat (rib) as followin 

i ribW 



(n b ) = 2 J d ® 



v /2(V($)-V(0) + e)/C 
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Supporting Material 



SI. Ideal crosslinking geometry of parallel actin bundles 

Here, we describe the details of the optimal packing calculation of our overtwist ground- 
state as shown in Fig. 2. The ideal crosslinking configurations should fulfill the following 
conditions: 1) maximize crosslinking bonds that are perfectly aligned to the hexagonal pack- 
ing directions and 2) minimize the twist distortion from the native —13/6 symmetry. To 
search for the optimal crosslinking configuration of parallel actin bundles, we first consider 
a single filament configuration with the commensurate helical geometry of —Qm/n sym- 
metry, in which each successive monomer is rotated by — 2irn/(6m), where both n and m 
are integers. This class of helical symmetry allows a subset of monomers align up with 
one of the sixfold lattice directions. In the block of m consecutive monomers with — 6m /n 
symmetry, two monomers at the boundary of the block make perfect alignments with the 
packing direction: for example, if the first monomer starts at <fi = 0, the mth monomer 
lines up with <fi = — 27rn/6. Among these helical symmetries, we consider the configurations 
parameterized by a pair of integers, (n, m), which minimize the angular deviation from the 
native symmetry. In Fig. ISTj we plot the angle difference from —13/6 symmetry for —Qm/n 
symmetry up to m = 7. As shown in Fig. IST1 it is found that —24/11 and —30/14 - re- 
spectively, under- and over-twisted relative to the native state - symmetry are closest to the 
native state, corresponding to m = 4 and m = 5, respectively. In our analysis, we focus on 
these two commensurate helical geometries. 

We now construct composite structures of a single filament, which consists of A/4 num- 
bers of 4-monomer sections with —24/11 symmetry and N 5 numbers of 5-monomer sections 
with —30/14 symmetry. For —24/11 symmetry, the 4-monomer section makes — 27r(ll/6) 
rotation, while for —30/14 symmetry, the 5-monomer section makes — 27r(14/6) rotation. In 
order to make full N turns, A/4 and N 5 must satisfy the following condition: — 4A r 4 (ll/24) — 
5AT 5 (14/30) = A" per M monomers, where M = AN^ + 5A5 and A^ is a negative integer for 
left-handed turns. As an example, the composite structures of a single filament with 28-fold 
symmetry are displayed in Fig. IS2j showing the longitudinal arrangements of blocks, such as 
445555 and 554455 in (a) and (b), and their top views with the orientations of crosslinkers 
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in (c) and (d), together with the crosslinker layer numbers. In table I, we present the class 
of composite structures, up to 102 monomers per repeat unit, that are used in this analysis. 

The remaining task is to build three-dimensional bundles by tiling the hexagonal lattice 
with the single-filament composite motifs for a given repeat length, M, and search for a 
state that maximizes the crosslinker density (defined as the ratio between the number of 
co-oriented monomers and the number of monomers in a bundle). Note that crosslinking 
only occurs between co-oriented monomers. For this task, we use a Monte Carlo method and 
find an upper bound of crosslinker density for a given symmetry. The initial configuration 
of the system is set by arranging filaments in the hexagonal array (5x5 rows). Each 
filament has one composite state of longitudinal block arrangements and one orientational 
state, which are randomly chosen from the block permutations of (iV 4 + N 5 )\ / N^. / N 5 \ states 
(all possible block translations and rearrangements along filament) and the six orientational 
states of hexagonal lattice directions. In the MC step, we make discrete rotational trial 
moves of filament by ±60° around its well as block translational trial moves, which 

are accepted if the trial moves increase the number of crosslinkers (i.e., the co-oriented 
monomers) and otherwise, rejected. The system is let to equilibrate until the crosslinker 
density is saturated, which usually requires from 1000 to 5 x 10 6 MC trial movements, 
depending on the number of monomers in the filament. To calculate an upper bound on the 
maximum crosslinker density allowed by a given composite structure, we count the number 
of all co-oriented monomers belonging to the 3x3 inner filaments, embedded within 5x5 
filament lattice. For a given finite size of lattice, this counting provides the upper bound 
of crosslinker density for an arbitrary large size of lattice, because we assume that this 
local packing arrangement may be continued over a larger region of the bundle. In fact, 
the crosslinkers along the boundary may be frustrated, resulting in an ultimately lower 
density of crosslinks upon considering a larger region of the bundle. Therefore, the number 
of co-oriented monomers for inner filaments serves as an upper bound for the maximum 
crosslinkers density allowed for a given symmetry. The MC procedure is repeated 1000 
times for a given N& and iV 5 to find the crosslinker configuration that gives the maximizes 
the number of co-oriented monomers. 

The upper bounds for crosslinker density are displayed in the chart of Fig. [S3] Up to 
102 monomer repeat units, we find that the groundstate geometry of 28-fold symmetry 
provides the maximum crosslinker density, that is, it saturates its upper bound with 6 
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bonds for every 28 monomers on every filament. However, for filament lattices of only 5x5 
rows, the upper bound on the crosslinker density of 13-fold symmetry is very close to that of 
optimal packing geometry. For the comparison with the optimal packing geometry of —28/13 
symmetry (see Fig. 2 in the text or Fig. IS6[) . we also display the crosslinking configuration for 

— 13/6 symmetry in Fig. ISB"l which has an upper bound of crosslinking density of 0.213675. 
Hence we performed more exhaustive searches on these two symmetries, by increasing the 
bundle size up to 10 x 10 to obtain tighter upper bounds for these structures, as shown 
in Fig. IS4l As the number of rows increases, the configuration with —28/13 composite 
filaments exhibits a constant crosslinkers density of 0.214286, while the crosslinking of the 

— 13/6 becomes increasingly frustrated, and ends up with a significantly reduced upper 
bound on the crosslinker density, compared to that of —28/13 structure, for 10 x 10 filament 
rows, the upper bound of the density of bonds is ultimately reduced to 0.198317 (see also 
the dark bar in Fig. |S3|) . The reduced value in the larger —13/6 bundles can be explained 
by the fact that, as the lattice size increases, the number of the "assumed" crosslinkers 
on boundary filaments decreases, which in turn reduces the upper bound. In the optimal 
crossliking geometry of —28/13 symmetry, unlike to the configuration with —13/6 symmetry, 
all the monomers aligned with the six-fold lattice directions are being fully consumed and 
involved in crosslinking. Note that —28/13 structure can tile a hexagonal lattice of arbitrary 
size, maintaining this maximum bond density. 

S2. 27r/28 rotational symmetry of the groundstate 

The optimal packing geometry of crosslinkning with —28/13 symmetry possesses the 
27r/28 rotational symmetry. Here, we demonstrate that under a coherent rotation by 2n/28, 
the bundle recovers the same crosslinking pattern of —28/13 overtwist groundstate, and 
hence, the free energy is degenerate under this rotation. Starting with the overtwist ground- 
state as shown in the left panel in Fig. IB6"| we take the following steps: 1) remove all 
crosslinkers, incurring a bond free-energy penalty, +AFb, 2) rearrange all the monomers 
into the homogeneous twist state of —28/13 symmetry by lowering twist energy, —AE tw . 
Notice that at this stage, monomers at layer and 14 are co-oriented along the horizontal 
lattice directions in Fig. [S6] Now, 3) rotate the entire filaments by 27r/28, which brings 
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monomers at layer 1 and 15 into perfect registry on neighbor filaments; 4) rearrange all 
monomers back to the original composite configurations by paying the same amount of 
twist energy that we gained, AE tw ; 5) replace the crosslinkers, regaining the same bond 
energy, —AF^. As a result of this coherent rotation, the layers associated with crosslinkings 
are changed as illustrated in the right panel in Fig. IB6l but the original crosslinking geometry 
and the total energy within 28 monomer repeat unit are preserved, albeit discreted rotated. 
Note that this entire procedure is identical to a one-layer translation along the bundle axis, 
followed by shifting the planar unit cell to the right by one-filament spacing. 



S3. Thermal fluctuations at finite temperature 



Our model of the commensurate-incommensurate transition described in the main text 
considers purely the elastic energy of filaments and crosslinkers, therefore it does not include 
the effect thermal fluctuations of monomers and crosslinker distributions. We here include 
these effects by considering all possible crosslinking geometries with a 28-monomer coarse- 
graining unit, and computing the probability of each state proportional to e~^ v ^ at finite 
temperature, /3 = ksT . We derive the effective pinning potential per monomer at the given 
temperature by considering the crosslinking statistics of each filament in bundle cross-section 
independently: 



Vf,($) = -k B T\n 



l + zJ2 exp(-/3K: ($)) + z 2 Yl e M-PV S2 ($)) 

si s 2 

+^^exp(-/3K 3 ($)) + ^ 4 ^exp(-/3K 4 (<f)) 

S 3 S 4 

+z 5 *M-PV S5 ($)) + z 6 *M-Pv a6 ($)) 



(SI) 

S5 S6 

where the fugacity is given by z = exp(/3(/i + e&)/2) and the summation runs over s nb , all 
possible bond distributions for given nj. Here, the potential for each state is defined by 



Vs. 



n b fi 



(S2) 



2 2 

Notice that the geometry of crosslinking requires us to consider a particular set of 19 corre- 
lated configurations of monomer orientations and crosslinker distributions in the 28-monomer 
span. In table II, we present the possible crosslinking configurations for given rib and the 
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corresponding the elastic offset energies, fi c (rib,s nb ). In Fig. [S7l we plot eq. (ISTj) revealing 
only modest quantitive differences in comparison to the "zero temperature" potential. These 
curves are quantiatively similar, though the sharp features associated with the changes in 
the minimal free energy bond pattern have been smoothed out by thermal fluctuations be- 
tween competing, nearly degenerate, bond configurations. To calculate the twist-state phase 
diagram (Fig. |S8|) including these finite temperature fluctuations, we set (3C = 6900 based 
on the known torsional stiffness of actin filaments C ~ 8 x 10 _26 Nm 2 per the monomer 
spacing of 2.8 nm. In Fig. (S9J the mean number of crosslinkers per 28-monomer span is also 
calculated by introducing the bond crosslinker density operator as follows: 

M ») = B6^ffi, (S3) 

along with the mean rotation profile, formula of eq. (B2). Including the thermal 

fluctuations of the distribution of bonds produces a phase-diagram that is quantitatively 
similar to the "zero-temperature" results presented in the text. Note, however, that the 
sharp features, kinks, in the (A0) and (rib) dependence that derive from abrupt changes in 
minimal-energy number state of crosslinkers that we saw in "zero temerapature" calculation 
(see Figs. 4 and 6 in the text) are smoothed out. In contrast, the sharp commensurate- 
incommensurate transition from the native to overtwist state is preserved at finite temper- 
ature, as the twist-suspectibility, d(A<f>) j dfi, necessarily diverges as the "lock-in" state is 
approached, indicating a second-order phase transition. 
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repeat length M 


helical symmetry 


overtwist (radians) 


(2,1) 


13 


-13/6 


0.000000 


(6,0) 


24 


-24/11 


-0.020138 


(2,4) 


28 


-28/13 


0.017261 


(4,5) 


41 


-41/19 


0.011788 


(10,2) 


50 


-50/23 


-0.009666 


(6,6) 


54 


-54/25 


0.008950 


(14,1) 


61 


-61/28 


-0.015847 


(12,3) 


63 


-63/29 


-0.007672 


(8,7) 


67 


-67/31 


0.007214 


(6,9) 


69 


-69/32 


0.014009 


(14,4) 


76 


-76/35 


-0.006359 


(10,8) 


80 


-80/37 


0.006042 


(20,1) 


85 


-85/39 


-0.01705 


(18,3) 


87 


-87/40 


-0.011111 


(16,5) 


89 


-89/41 


-0.005431 


(12,9) 


93 


-93/43 


0.005197 


(10,11) 


95 


-95/44 


0.010175 


(8, 13) 


97 


-97/45 


0.014948 


(22,2) 


98 


-98/45 


-0.014796 


(18,6) 


102 


-102/47 


-0.004738 



TABLE I. The composite structures of a single filament constructed by and N§ blocks with 
—24/11 and —30/14 symmetry, respectively. We analyze up to 102 monomer repeat length. 
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FIG. SI. The angle difference per monomer of —6m/n symmetry from the native —13/6 symmetry. 
— 24/ll(m = 4) and — 30/14(m = 5) symmetry are nearest to the native state of twist. The red 
dashed line is the angle deviation of —28/13 symmetry. 
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FIG. S2. The structural details of composite structure for 28-monomer repeat with N4 = 2 and 
N$ = 4. The side views of filament composed by block arrangements of 445555 and 554455 are 
shown in (a) and (b), respectively and their top views are shown in (c) and (d). The layer numbers 
are the location of crosslinkers. 
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FIG. S3. The upper bound of crosslinker density for various symmetries up to 102 monomers 



are presented in bar charts. The results are searched from a Monte Carlo calculation for 5x5 



bundle size. The 28-fold symmetry is found to be the structure with maximum crosslinker density, 
6/28. For 13-fold symmetry, the upper bound crosslinking density, 0.198317, from 10 x 10 row MC 
calculation is marked with a dark gray bar. 
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FIG. S4. The explicit comparison between —13/6 and —28/13 symmetry with respect to the row 
number in MC bond calculation. Note that 28-fold symmetry shows a constant crosslinker density 
(saturating the upper bound, 6/28), while the upper bound decreases for the 13-fold symmetry as 
the row number grows. 
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FIG. S5. The maximum crosslinking configuration for —13/6 symmetry by Monte Carlo calculation. 



We calculate an upper bound on the crosslinker density by counting crosslinks formed inside the 



dashed line. Note that the crosslinkers at the boundary are share with filaments outside the cell, 



contributing only half a crosslink to the upper bound counter. For a 5 x 5 filament structure, the 



upper bound of crosslinker density is (6 x 2 + 13)/(13 x 9) = 0.213675. 
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5,18 4,15 -57T8-- 4,19 5,20 fi"l9 5 20 " 6 . 19 

FIG. S6. The crosslinking structure in the right is generated by the 2-7r/28 coherent rotation 
of maximally crosslinked bundle groundstate shown in the left, followed by a reorganization of 
monomers and crosslinks within a 28-monomer repeat length. Note that under rotation the layers 
where crosslinkings occur shift, while the crosslinking directions are preserved. 
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n b = 


so 


1 28 | 












-e& 










n b = 1 


Si 


1 28 | 












-e& 










n fe = 2 


S2 


| 14 | 14 | 


| 5 | 23 | 


1 8 | 20 | 


| 10 | 18 | 


| 13 | 15 | 




-e& 


-e b + 0.00068C 


-e b + 0.00783C 


-e b + 0.00174C 


-e b + 0.00361C 


ra 6 = 3 


«3 


| 4 | 10 | 14 | 


| 9 | 5 | 14 | 


| 13 | 5 | 10 | 


| 8 | 5 | 15 | 






-e 6 + 0.00261C 


-e b + 0.00058C 


-e b + 0.00241C 


-e b + 0.00522C 




n b = 4 


S 4 


| 4 | 5 | 5 | 14 | 


| 4 | 10 | 5 | 9 | 


| 4 | 10 | 10 | 4 | 


| 9 | 5 | 5 | 9 | 


| 8 | 5 | 10 | 5 | 




-e 6 + 0.00195C 


-e b + 0.00239C 


-e b + 0.00391C 


-e b + 0.00087C 


-e b + 0.00391C 


ra fe = 5 


S5 


| 4 | 5 | 5 | 5 | 9 | 


| 4 | 5 | 5 | 10 | 4 | 










-e b + 0.00191C 


-e b + 0.00313C 








n fe = 6 


S6 


|4|5|5|5|5|4| 












-e b + 0.00261C 











TABLE II. The crosslinking configurations commensurate with the -28/13 filament geometry for 
all n b and the corresponding elastic offset energy, fi c (n b , s nb ). 
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FIG. S7. The $ dependence of the crosslinking free energy at f3C = 6900 for U = 3.0C (top) and 
U = 0.5C (bottom) for various chemical potentials. Finite temperature effects from the statistical 
distribution of crosslinks in the bundle lead to modest quantitative differences with Fig. 3 of the 
text. 
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FIG. S8. The equilibrium diagram of state for overtwist, (A0), for (3C = 6900. The CI transition 
is shown as a solid line. 
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FIG. S9. The mean number of crosslinker per 28-monomer repeat, for f3C = 6900. The solid 
line denotes the CI transtion. 
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